Sulphur abundances from the Si near-infrared triplet at 1045nm ^ 
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ABSTRACT 



Context. Unlike silicon and calcium, sulphur is an cy-element which does not form dust. Some of the available observations of the evolution 
of sulphur with metallicity indicate an increased scatter of sulphur to iron ratios at low metallicities or even a bimodal distribution, with some 
stars showing constant S/Fe at all metallicities and others showing an increasing S/Fe ratio with decreasing metallicity. In metal-poor stars S i 
lines of Multiplet 1 at 920 nm are not yet too weak to permit to measure the sulphur abundance A(S), however in ground-based observations 
they are severely affected by telluric lines. 

Aims. We investigate the possibility to measure sulphur abundances from S i Mult. 3 at 1045 nm lines, which lie in the near infra-red. These are 
slightly weaker than those of Mult. 1, but lie in a range not affected by telluric lines. 

Methods. We investigate the lines of Mult. 3 in the Sun (G2V), Procyon (F5V), HD 33256 (F5V), HD 25069 (G9V) and 6 Eri (HD 22049, 
K2V). For the Sun and Procyon the analysis has been performed with CO^^BOLD 3D hydrodynamical model atmospheres, for the other three 
stars, for which hydrodynamical simulations are not available, the analysis has been performed using ID model atmospheres. 
Results. For our sample of stars we find a global agreement of A(S) from lines of different multiplets. 

Conclusions. Our results suggest that the infrared lines of Mult. 3 are a viable indicator of the sulphur abundance which, because of the intrinsic 
strength of this multiplet, should be suitable to study the trend of [S/Fe] at low metallicities. 
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1. Introduction 

The so-called ff-elements (O, Ne, Mg, Si , S, Ca) are among 
the main products of T ype II Su pernovae (^Woosley & Weavei 
ll995l:lLimongi&Chieffi 2003a; Chieffi&Limongi 2004). The 
iron peak elements are produced by Type II SNe, and it is com- 
monly accepted that Type la Supernovae produce similar, or 
larger amou nts of iron-peak elements, and produce li ttle or no 
a-elements, dNomoto et all 1 984tllwamoto et al J 1 999i) . The dif- 
ferent timescales for the first explosions of Type II or Type la 
SNe to occur, makes the abundance ratio of c-elements to iron- 
peak elements a powerful diagnostics of the chemical evolution 
and star formation history of a galaxy. In the Milky Way, stars 
of lower metallicity are characterised by higher a-to-iron abun- 
dance ratios than found in the Sun and solar metallicity stars 
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(see, for example, Barbuy 1988; Gratt on et al.l 120031; Cavrel 
et aI. T2004,) . This is usually interpreted in terms of a lower con- 
tribution of Type la SNe. Systems which are characterised by 
low or bursting star formation, like dwarf galaxies, give time to 
Type la SNe to explode before the enrichment due to Type II 
SNe has greatly increased. Consequently, such systems display 
rather lo w g-to-iron ratios even at low metallicities (se e Venn 
et al. i200i" and referen ces therein) and at solar metallicities . 



displa y sub-solar ratios dBonifacio et al.ll2004t [Monaco et al 



2005) 



One should be aware that the above outlined, simple inter- 
pretation of the a-to-iron peak element ratios in terms of prod- 
ucts of Type II and Type la SNe, relies on the nucleosynthesis 
computed with ID explosion models. For Type la SNe such 
computations predict that the original C-0 white dwarf is to- 
tally burned, mainly to ^^Ni . It is i nteresting to note that in the 
2D models by Brown et alj d2005h . most of the white dwarf is 
not burnt and roughly equal masses of ^^Ni and ^^Si are pro- 



duced. It IS clear that a dirrerent explanation or the c-element 
enhancement must be sought if these computations are con- 
firmed. 

While the theoretical interpretation of the a-to-iron peak 
ratios is likely still open to debate, it is clear that, from the ob- 
servational point of view, the abundance of a-elements is an 
important property of any stellar population. For the study of 
chemical evolution in external galaxies, the more readily avail- 
able objects are Blue Compact Galaxies (BCGs) through anal- 
ysis of their emission line spectra, and Damped Ly-a systems 
(DLAs) through the analysis of resonance absorption lines. In 
this way, it is relatively easy to measu re sulphur in t he gaseous 
component o f both groups of galaxies dOarnettl 19891: Centurion 
et al. ,2000) . 

At variance with Si and Ca, S is a volatile element, and 
therefore it is not locked into dust grains in the interstellar 
medium, so that no correction is needed to the measured sul- 
phur abundance. This makes sulphur a more convenient ele- 
ment to trace the a's than either Si or Ca; all three elements 
are made by oxygen burning, either in a central burning phase, 
conve ctive shell, or explosive phase. According to Limongi & 
Chieffi (1200^ there is thus a strong reason to believe that 
Si, S, and Ca vary in lockstep during the chemical evolution. 
In spite of this, it is certainly unsatisfactory to compare S/Fe 
in external galaxies to Si/Fe or Ca/Fe in the MiUcy Way. It is 
much more convenient to use a reliable Galactic reference for 
sulphur abundances, which may be directly compared to mea- 
sures in external galaxies. 

The only way to measure abundances at different metal- 
licities in the Galaxy is to use stars. Unfortunately there are 
very few sulphur lines which are not blended and remain 
strong enough to be measured at low metallicities. The lines 
of Mult. ^ (675 nm) and 6 (869 nm) are weak, so only de- 
tectable in solar or moderately metal-poor stars, hardly lower 
than [Fe/H] ~ -1.5 (for the Muh. 8), or -2.0 (for Mult. 6). The 
lines of Mult. 1 (920 nm) have been recently used to measure 
A(sil because these lines are strong, and are detectable even 
at low metallicity. The non-ne gligible NLTE effec ts in these 
lines make them even stronger jTakeda et al ] l2005ah . therefore 
more easily measurable. The main problem is that the range in 
wavelength in which Mult. 1 lies is contaminated by numerous 
telluric lines. This makes it difficult to obtain all (or at least 
one) components of Mult. 1 unaffected by telluric absorption. 

The 1045 nm lines of Mult. 3 are particularly suited to mea- 
sure the sulphur abundance. Even if the lines are not as strong 
as the components of Mult. 1, the big advantage is that there 
are no telluric lines present in the vicinity of their wavelength. 
Observing sulphur lines of Mult. 3 provides a possibility to ob- 
tain a reliable sulphur abundance in very metal-poor stars. 

Existing studies of sulphur in Galactic stars have suggested 
that in the range -2.5 <[Fe/H]< -2.0 the [S/Fe] ratio shows 



either a large scatter or a bimodal behaviour; Most of the stars 
converge to a "plateau" at about [S/Fe] =-1-0.4, while a non- 
negligible number of stars s hows a "high" va lue of [S/Fe] 
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around -1-0.8. In the sample of Caffau et alJ (120051) . the determi- 
nation of [S/Fe] in this range of metalUcity is based on the non 
contaminated lines of Mult. 1 . This behaviour has no proposed 
theoretical interpretation and it may well be due to systematic 
errors related to the use of Mult. 1 . It is therefore of great inter- 
est to verify this puzzling finding by the use of an independent 
and, hopefully better, diagnostics of the sulphur abundance, as 
can be afforded by the lines of Mult. 3. We moreover recall that 
these are the only strong S i lines belonging to a triplet system 
instead of quintet systems as the other S i lines in the visual and 
near-IR range. 

The aim of the paper is to study the S i lines of Mult. 3 in 
two well known stars - the Sun and Procyon - and to investi- 
gate three other bright field stars of spectral type F, G and K: 
HD 22049, HD 33256, HD 25069 respectively. We compare 
A(S) from two sulphur lines widely used in sulphur abundance 
determinations, to the one from Mult. 3 lines. We want to es- 
tablish Mult. 3 as valuable abundance indicator. 



2. Atomic data 

The S I lines that we consider in this work, are reported in Table 
[U The triplet at 675.7 nm and the line at 869.4 nm have been 
widely used in the determination of sulphur abund ances (see 
for instance. iNissen et al.ll2004l ; ICaffau et al.ll2005h . The lines 
of Mult. 3 have not yet been consider ec E ex cept in the Sun 
jLambert & Lucklll978l: iTakeda et al.ll2005d) . We know that 
the bluest line of Mult. 3 is blended with an iron line of poorly 
known loggf, but we nevertheless keep this line for the Sun 
and Procyon for which we have good observed spectra in hand, 
and for HD 33256 where A(Fe) is lower than in the Sun, as well 
as for e Eri where the iron contribution is only 12 % of the total 
equivalent width (EW); we discard this line for HD 25069 for 
which the contribution of iron to the total EW is about 18 %, 
comparable to the error of the EW for this star In any case, 
the relative contribution of the iron line to the blend becomes 
smaller for very metal-poor stars, and thus the line can be a 
good indicator of the sulphur abundance. 

The logg/ value of the Mult. 3 lines, as those of 
all the other S i lines used here, have been taken from 
the Kurucz line list; the data are approximately coinci- 
dent with those from the NIST Atomic Spectra Database 
(|http://physics.nist.gov/PhysRefData/ ADS[ ). These gf are ex- 
perimental values based on three experi ments. Owing to the 
moderate accuracy of the absolute scale, IWiese et alJ (Il969h 



We use the Multiplet numbering bv lMoorel ( [T94^ 

2 A(S)=log (N(S)/N(H)) + 12 



The referee made us aware of a paper (Nissen et al. 2007 astro- 
ph/0702689, submitted to Ac&A) which became available through 
arXiv (arxiv.org) after the submission of this work. The paper de- 
scribes the first use of the 1045 nm lines for abundance work in metal- 
poor stars. 



lable 1. Atomic parameters or the sulphur hnes. Col. (1) is the 
wavelength; col. (2) the multiplet number; col. (3) the transi- 
tion; col. (4) the log g f of the transition taken from the Kurucz 
line list; col. (5) the excitation energy; col. (6) the Van der 
Waals damping constant computed at a temperature of 5 500 K 
according to ABO theory (see Sect. |6]l or Kurucz approxima- 
tion (see text). 



Wavelength 
(nm) air 
(1) 



Mult. Transition log gf 



(2) 



(3) 



(4) 



(eV) 
(5) 



logye/Nn 
(6) 



675.6851 
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-1.76 


7.87 


-7.146 


675.7007 
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-0.90 


7.87 


-7.146 


675.7171 
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-0.31 


7.87 


-7.146 


869.3931 


6 






-0.51 


7.87 


-7.337 


869.4626 
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^^P3 




0.08 


7.87 


-7.337 


1045.5449 


3 


3S° 


-3p, 


0.26 


6.86 


-7.672 


1045.6757 


3 




-^Po 


-0.43 


6.86 


-7.672 


1045.9406 


3 




-^Pi 


0.04 


6.86 


-7.672 



renormalised the data theoretically and judged the resulting ac- 
curacy of the gf to be within 50 %. In particular we note that all 
the lines selected for S abundance determination are of D qual- 
ity, i.e. the uncertainty of the oscillator strength is < 50 %; the 
lines with highest accuracy (D-n, i.e. an uncertainty of the os- 
cillator strength < 40 %) are those of Mult. 1 (920 nm), Mult. 3 
(1045nm), and the Hne at 869.4nm of Mult. 6, while the Hne 
at 869.3 nm of Mult. 6 is of D quality, as is one component 
(675 .7 171 nm) of the 675 .7 nm line of Mult. 8 . One component 
of the 675 nm triplet of Mult. 8, the 675.6851 nm Hne, is of E 
quality (uncertainty of the oscillator strength > 50 %). 

3. Models 

For the Sun and Procyon our analysis is based on 3D hydro- 
dyna mical model atmosphe res computed with th e CO^BOLD 
code JFreytag et al . 2002; Wedemeyer e"taDl2004l) . CO^BOLD 
solves the coupled non-linear equations of compressible hy- 
drodynamics in an external gravity field together with non- 
local frequency-dependent radiation transport for a small vol- 
ume located at the stellar surface (see CO^BOLD manual 
http://www.astro.uu.se/~bf/cobold/index.html). 25 snapshots 
were selected from a CO^BOLD simulation to represent the 
photosphere of the Sun (Caffau et al., accepted for A&A 
Letters), with an effective temperature of 5780 K and cover- 
ing 6000 s of temporal evolution. 28 snapshots were selected 
from a 3D simulation of Procyon, having an effective tem- 
perature of 6500 K and covering a time interval of 16800 s. 
Since the timescale of the evolution of the granular flow is 
about 2.8 times longer in Procyon than in the Sun, the sim- 
ulated time spans are very similar in a dynamical sense. The 
Procyon model used in this paper is the same as the one used 



bv lAutdenberg et al.l 020051) in their study of the star s centre- 
to-limb variation. In both hydrodynamical models, an opac- 
ity binning scheme with five wavelength bins was applied for 
modelling the wavelength-dependence of the radiative transfer 
(■Nordlund.l982 ;.Ludwig et al..l994; Yogler et al...2004) . 

For the other stars we do not have available 3D atmo- 
spheres, and therefore used ID LTE plane-parallel models. 
Such ID models were also used as a reference for the Sun and 
Procyon. In particular: 

1. For all the stars we computed hydrostatic ID model 
atmospheres computed with the LHD code. LHD is a 
Lagrangian ID (assuming plane-parallel geometry) hydro- 
dynamical model atmosphere code. It employs the same 
micro-physics (equation-of-state, opacities) as CO^BOLD. 
The convective energy transport is described by mixing- 
length theory. The spatial discretisation and numerical so- 
lution of the radiative transfer equation is similar to the one 
in CO^BOLD, albeit simplified for the ID geometry. The 
wavelength-dependence of the radiation field is treated by 
the opacity binning method. A hydrostatic stratification in 
radiative-convective equilibrium is obtained by following 
the actual thermal and dynamical evolution in time of the 
atmosphere until a stationary state is reached. LHD pro- 
duces standard ID model atmospheres which are differen- 
tially comparable to corresponding 3D CO^BOLD models. 
Remaining choices entering an LHD model calculation are 
the value of the mixing-length parameter, which formula- 
tion of mixing-length theory to use, and in which way tur- 
bulent pressure is treated in the momentum equation. Note, 
that these degrees of freedom are also present in other ID 
model atmosphere codes. The LHD models presented in 
this work have all been computed with ctmlt = 15 using 
the formulation of iMihalasI (Il978h . and turbulent pressure 
has been neglected. Comparisons to 3D models were al- 
ways made with LHD models having the same effective 
temperature, gravity, and chemical composition as the 3D 
model. 

2. For all stars, except the Sun, we computed ATLAS9 
Kuruczlll993aL al) models using t he Linux version 



Sbordone et al. 



20041; ISbordond 120051) of the code. AU 



these models have been computed with the "NEW " Opacity 
Distribution Functions ( Castelh & Kuruczl2003h. which are 



Grevesse & Sauvalld 1998b 



based on solar abundances from 

' micro-turbulence, a mixing-length parame 



with 1 km s 
ter qtmlt of 1 .25 and no o vershooting. 



3. We used an ATLAS9 (iKuruczl Il993al l2005al) model 
of the Sun comput ed by Fiorella Castell i with the so- 
lar abundances of lAsplund et al.l (l2005h as input for 
the chemical composition. The Opacity Distribution 
Functions were explicitly computed for the same 
chemical composition and a micro-turbulent velocity 
of Ikms"'. The model was computed assuming a 
mixing-length parameter of 1.25 and no overshooting. 



It IS available at http://wwwuser.oats.inar.it/castelli/sun/ 
ap00t5777g44377klasp.dat. 



We used the Holweger-Miiller solar model (iHolweger 
1967l:lHolweger & MuelleJl974 . It is an empirical model 



of the solar photosphere and lower chromosphere assum- 
ing LTE, and largely reproducing the Sun's continuous and 
line spectrum. It considers 900 selected line profiles of 3 1 
atoms and ions. 

5. For the Sun and Procyon we considered horizontal and tem- 
poral averages of the 3D snapshots over surfaces of equal 
(Rosseland) optical depth. Comparison with these averaged 
3D models, henceforth denoted as (3D) models, provide 
estimates of the influence of fluctuations around the mean 
stratification on the line formation process. Comparing 3D 
and ID models of this kind is largely independent of arbi- 
trary assumptions entering the constructions of standard ID 
models. The only free parameter which has to be specified 
for the ID average model is the micro-turbulence to be used 
in related spectrum synthesis calculations. 



ettects, as would be the case it we computed corrections with 
respect to other hydrostatic models, like ATLAS. It may be use- 
ful to summarise the differences between a 3D hydrodynamic 
model and a ID hydrostatic model as due to two factors: 

1. a different mean structure of the two classes of models; 

2. horizontal temperature and pressure fluctuations, which are 
present in a 3D simulation but not (by definition) in a ID 
model. 

Our definition of "3D correction" is made in order to take into 
account both effects at the same time. 

If one wants to see the difference in abundance due only to 
the different mean temperature structure one may compare the 
abundances derived from the horizontally and temporally aver- 
aged 3D model, which we defined above as the (3D) model, to 
the corresponding LHD model. On the other hand, if one wants 
to evaluate only the effect of horizontal temperature and fluc- 
tuations, it is more appropriate to compare the 3D abundances 
with those derived from the (3D) model. 



The spectral synthesis calculations for the 5_ Qata 



CO^BOLD, LHD, ATLAS and Holweger-Muller mod- 
els were performed with the code Linfor3D (see 



http://www.aip.de/~mst/Linfor3D/linforJ3D_manual.pdf); 



For the Sun we considered two high resolution, high signal to 
noise ratio, spectra of the solar flux. 



for AT LAS models we also used SYNTH E ("Kurucz" 1993b', 
2005 al) . in its Linux version dSbordone et al. 2004 ; Sbordone 
20051) . This was used in the cases in which we wanted to 
compute a synthetic spectrum containing many lines from 
other atoms and molecules, to fit the observed data. The 
present version of Linfor3D can handle at maximum a few 
tens of lines at a time, while SYNTHE does not suffer such a 
limitation. 

4. 3D abundance corrections 

The purpose of analysing a star both with 3D hydrodynami- 
cal models and ID models is to derive "3D abundance cor- 
rections" which may be used to correct the analysis of other 
stars of similar atmospheric parameters, performed in ID. It is 
clear that from the computational point of view ID modelling 
(both model atmosphere and spectrum synthesis) is much more 
convenient than 3D modelling. It is unlikely that in the near 
future the computing power will increase to the level that 3D 
analysis of stellar spectra will be done routinely. On the other 
hand it is conceivable that grids of hydrodynamical models will 
be constructed from which "3D corrections" can be computed, 
and these can be used to correct the results of a ID analysis. 
It is therefore important to define what is meant by a "3D cor- 
rection". In this paper we give 3D- ID abundance corrections 
with reference to LHD models. These provide a reasonably 
well-defined way to establish the relation between a 3D hy- 
drodynamical and ID hydrostatic model. The fact that LHD 
and CO^BOLD use the same opacities and micro-physics en- 
sures that the differences reflect only the 3D effects and no other 



1. 



The spectrum we will refer to as "Kurucz flux" is based 
on fifty solar Fourier transform spectrometer (FTS) scans 
taken by J. Brault and L. Testerman at Kitt Peak, with a 
spectral resolution of the order of 300000 and signal to 
noise of around 3 000, varying from range to range (further 

details can be found in Kurucz 200 5B. 

The "Neckel flux" refers to the Nec kel & Labsl (11984 1 abso- 



lutely calibrated FTS spectra obtained at Kitt Peak, cover- 
ing the range 330 nm to 1250 nm. The spectral purity ranges 
from 0.4pm at 330 nm to 2pm at 1250nm. This means that 
the resolution at 1045 nm is about 500000. 

The spectra for Procyon and the other stars were obtained 



from the UVES Paranal Observatory Project (IBagnulo et al 



20031). For the 670 nm and the 870 nm lines, we have taken the 
reduced data present in the UVES POP web sitfl the spec- 
tral resolution is of about 80 000. The lines of Mult. 3 are not 
available in the POP reduced spectra, although inspection of 
the raw data reveals that these wavelengths are indeed recorded 
on the MIT CCD of UVES in the standard 860 nm setting, but 
the number of counts is very low (UVES is very inefficient at 
these wavelengths) and only about 1/3 of the order is present 
on the CCD. For these two reasons a standard extraction with 
the UVES pipeline fails to extract this order We downloaded 
the raw data and associated calibration observations from the 
ESO archive and reduced them using the UVES pipeline. We 
forced the extraction of the last order in the 860 nm setting by 
declaring the number of orders to be found. This allowed a sat- 
isfactory order definition and subsequent order extraction. 



http://www.sc.eso.org/santiago/uvespop7l 



b. uata analysis 

To derive sulphur abundances we measured the equivalent 
width (EW) of the selected sulphur lines and, when possible, 
fitted the line profiles. 

For the triplet lines, EWs were computed by direct numer- 
ical integration using the trapeze sum rule, or with the inte- 
gration of a fitted Gaussian for weak Unes or Voigt profile for 
strong lines, using the IRAF task splot. 

All line profile fitting has be en performed using a code. 



proiile iittmg 

described in ICaffau et al.l (120051) . which performs a rnin 



imisation of the deviation between synthetic profiles and the 
observed spectrum. In the fitting, the sulphur abundance, the 
level of the continuum, and a wavelength shift are left as free 
parameters, while the macro-turbulence has to be fixed a priori. 

When available, for the Van der Waals broadening we used 
the pa r ameter derived frorn the theory bvlAnstee & O'Mara 
(ll995h : lBarkle m & 0'Mara'('l99'/):'B arklem et al.l(ll998b l) and 
summarised in.Barklem et al. ( 1998a^ ); in this paper we collec- 
tively refer to this work as "ABO" theory. Otherwise we relied 
on the approximation build into the SYNTHE spectrum syn- 
thesis suite (Kurucz 1993 b, 2005a). The co mparison of t his ap- 
proximation to other ones can be fo und in iRvanI (1 19981) . who 
refers to its use in the WIDTH code ("Kurucz '2005a'). We note 
here that the approximation is used in SYNTHE and WIDTH 
for all lines for which literature data on the Van der Walls 
broadening do not exist, and not only for iron-peak elements, 
as the reader might be induced to believe reading the paper of 
Rvanl(ll998h . 

We used a 3D CO^BOLD simulation for the sulphur 
analysis in the Sun and Procyon. For e Eri (HD 22049) 
we computed the 3D- ID abundance corrections by using 
a CO^BOLD simulation whose parameters (reff=5073K, 
log ^=4.42, [M/H]=0.0) were very close to its stellar param- 
eters. 



6.1. Corrections for departure from local 
thermodynamic equilibrium 

According to iTakeda et al.] (l2005al) several of the atomic sul- 
phur lines show non-negligible departures from LTE. It is be- 
yond the purpose of this paper to investigate these departures, 
howev er we shall make use of the published departures from 
LTE of iTakeda et all (l2005ah . 



7. Results for the individual stars 

7.1. Sun 

For the Sun we consider the 675.7 nm triplet of Mult. 8, the 
869 nm lines of Mult. 6 and the three lines of Mult. 3 even if 
we know that the bluest is blended. 

The results are reported in Table |2] and are: the EW for 
each feature (second column), the sulphur abundance derived 
from the EWs and the results for the 3D corrections, assuming 



a micro-turbulence fmicro of 1 .5 km s" and 1 .0 km s" m the ID 
models. The 3D corrections are always defined with respect to 
the corresponding ID LHD model. In LTE, the 3D abundance 
corrections for the Sun turn out to be small. By magnitude this 
is compatible with the findings of Nissen et al.. QQQA), how- 
ever, our corrections have opposite sign. It is quite likely that 
the difference comes about by the different choice of ID refer- 
ence atmosphere to which the corrections are related. In fact. 



Nissen et all (120041) used a MARCS model as the ID refer- 



ence. In the table the NLTE correction, accordina to Takeda 
et al. (2005a), is included in the last but one column (A). The 
last column is the adopted sulphur abundance. In columns 4-11 
there are the sulphur abundance values derived from ID mod- 
els, using a^micio of LO and L5kms"'. 

In Fig. [T] the solar sulphur lines considered in this work 
are plotted together with the 3D synthetic spectra. The 
weak lines of Mult. 6 and Mult. 8 have also been fitted 
with both the CO^BOLDH-Linfor3D synthetic spectra and 
ATLAS9+SYNTHE. The results are, within 0.05 dex, in agree- 
ment with the abundances obtained from the EWs. 

The LTE computation implies A(S)=7.252 + 0.143 (this is 
an average of all the A(S) given in column (3) of Tabled, while 
applying the NLTE corrections of ITakeda et al.l (12005 a) (given 
in column (11) of TableEJ we obtain A(S)=7.213 ± 0.113; the 
quoted errors are the line to line scatter in the abundance deter- 
minations from the six features considered. It is worth pointing 
out that, given the high S/N of the available solar spectra, the 
associated statistical errors are negligible. Therefore the line 
to line scatter must reflect inadequacies in the model atmo- 
spheres and/or the line formation calculations, and/or errors in 
the atomic data. 

Examining more closely the disagreement between lines, 
we note that the 869.3 nm line is blended with molecular lines 
(C2 and CN) whose contribution we estimated to be of the or- 
der of 15 %, and which was consequently subtracted from the 
EW. However, considering the uncertainties in the logg/ of 
these molecular lines, it is probably safer to discard this line. 
From the 675.7 nm triplet we find A(S)=7.136 + 0.004 (where 
the error is now the standard deviation between the measures 
from the two observed solar spectra), while the sulphur abun- 
dance is 7.204 ± 0.021 in LTE and 7.194 + 0.021 in NLTE for 
the 869.4 nm fine. From Mult. 3 we find A(S)=7.379 ± 0.014 
in LTE and A(S)=7.309 ± 0.016 in NLTE. There is an obvious 
trend in A(S) with EW. Note that since this result is obtained 
using 3D atmospheres we cannot invoke a micro-turbulent ve- 
locity to remove this trend. This eff'ect could be explained if 
3D-NLTE corrections are larger than the published ID-NLTE 
corrections, adopted here. If this were the case, the solar S 
abundance should be A(S)=7.14, as indicated by the 675.7 nm 
triplet, which is virtually unaffected by NLTE. This result is 
also supported by the analysis of [SI] line at 1082 nm (Caff'au 
& Ludwig, in press), whose departure from LTE is negligible, 
from which the sulphur abundance is 7.14. However, the trend 
could have other explanations. At least in part, the line broad- 



lable 2. Solar sulphur abundances rrom tlux spectra. Col. (1) is the wavelength or the line rollowed by an identihcation Hag, 
K means Kurucz flux, N Neckel flux; col. (2) is the Equivalent Width; col. (3) is the sulphur abundance, A(S), according to the 
CO^BOLD 3D model; col. (4) to (11) give the abundances for a micro-turbulence of 1.5 and 1.0km s ' for the <3D), ATLAS, 
HM, and LHD model, respectively; col. (12) and (13) provide the 3D abundance corrections relative to the LHD model for a 
micro-turbulence of 1.5 and 1.0 kms"', respectively; col. (14) is the NLTE correction from .Takeda et al., (.2005a) ; col. (15) the 
A(S) we finally adopted. LHD models were computed with cmlt of 1.5, ATLAS models with aMLT of 1.25. 
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Fig. 1. Sun: synthetic spectra (green/grey thick solid line) based on the 3D model atmosphere are over-imposed on the observed 
(black solid line) solar flux spectrum of Kurucz. 
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Fig, 2. Procyon: the 3D fit (green/grey crosses) is over-imposed 
on the observed spectrum (black soUd line) for the S 1 675.7 nm 
triplet. 

ening theory employed could be inadequate, as could be the 
case for the hydrodynamic velocity and tempe rature field of 



the CO BOLD simulation. Finally, we note that lAsplund et al 



(I2005h find a solar LTE sulphur abundance of A(S)=7.14±0.05. 



7.2. Procyon 



For Pr ocyon we adopt the stellar parameters derived bv lSteffen 
(Il985h . that is T^g = 6500 K, logg=4.0 and solar metalhcity. 
These stellar parameters are very close to the results o f the most 
recent determination due to lAufdenberg et al. I (l2005h . who find 
reff=6516 + 87 K, log ^=3.95 + 0.02. For Mult. 8 and Mult. 6 
we measure S/N=400, for Mult. 3 S/N=50. 

For the 3D fitting we assumed a Gaussian instrumental 
broadening with FWHM of 3.75 km s ', corresponding to a 
spectral resolution of 80 000. 

The results are reported in Table [3] In the second column 
the EW is given, the next five columns show the sulphur abun- 
dance computed from EWs. The next two columns show the 
sulphur abundance from line profile fitting (see in Fig. |2]the 
fit of the 675. 7nm triplet); the two lines (869.3931 nm and 
869.4626 nm line) of Mult. 6 are fitted at the same time. In 
column (10) the 3D- ID correction is given, in column (11), 



when av ailable, the NLTE correction according to lTakeda et al 



(l2005ar) is reported. In the last but one column the adopted sul- 
phur abundance is given. 

The standard deviation of the EW relat ed to th e finite S/N 
is computed according to Cayrel's formula ( Cayrell 1988.) : 



e(EW) = 1.6 X 



VFWHM X PixelSize 



S/N 



(1) 



where S/N is the signal to noise value, FWHM is the full width 
at half maximum of the line, PixelSize is the size of the detec- 
tor pixel in wavelength units. From the errors of the EWs we 
computed the corresponding errors of A(S) using the curve of 



growth of each Ime. Ihese errors are provided m the last col- 
umn of Table [3] 

For Procyon the 3D corrections are definitely not negligi- 
ble, above all for the strong lines of Mult. 3. This effect is not 
so evident in the Sun which is cooler than Procyon. We shall 
see that the 3D corrections are negligible in the cooler e Eri 
(see Table |6]l. We note that the 3D corrections given in Tab. [3] 
are positive, implying that - for a given abundance - the same 
fines are weaker in 3D than in ID. This beh aviour is opposite 
to the findings by lSteffen & Holwegen(l2002h who argue that in 
general one should expect a strengthening of lines in 3D atmo- 
spheres. However, their argumentation is not strict, and actually 
refers to the 3D-(3D) correction for weak (unsaturated) lines in 
the Sun. In fact, the weaker lines in Tab. |2] comply with the 
expected line strengthening in 3D. 

The fact that the difference of the 3D abundance to that 
derived from the (3D) model (col. (3) - col. (4) in Tab. O is 
of similar size as the total 3D- ID correction given in col. (9) 
(col. (3) - col. (6) in Tab. |3]l shows that the pronounced line 
weakening cannot be the result of a different mean temperature 
structure in 3D with respect to ID. In fact the abundances de- 
rived from the (3D) (col. (4)) and the LHD model (col. (6)) are 
very similar Rather the line weakening must be caused by the 
sizable horizontal temperature fluctuations present in the pho- 
tosphere of Procyon, or by a substantially different effective 
micro-turbulence in 3D and ID, respectively. 

In order to explore this question somewhat further, 
Figures|3]and|4]depict the joint probability density of the (disk- 
centre) EW of the 1045.9 nm line and the neighbouring contin- 
uum intensity for the Procyon and solar 3D model, respectively. 
The contours illustrate the correlation between continuum in- 
tensity and EW over the stellar surface. We find a bimodal dis- 
tribution in the solar case, a single peak and diffuse "halo" in 
the case of Procyon. The difference can be traced back to quali- 
tatively different formation heights of the line. In the case of the 
Sun the largest contribution to the line absorption stems from 
the granular layers as such, while the maximum contribution 
is shifted to higher layers of reverse granulation in Procyon. 
While in the Sun temperature fluctuations in the continuum 
formation layers are positively correlated with the tempera- 
ture fluctuations in the line forming layers, we typically find 
an anti-correlation in Procyon. Due to the high excitation po- 
tential of the lower level the line is rather temperature sensitive. 
In the Sun instances of high continuum intensity coincide with 
a high temperature in the line forming layers. Higher tempera- 
tures increase the population of the lower level of the transition 
leading to a stronger line. The reverse happens in moments of 
low continuum intensity. In Procyon, the line still tends to be- 
come stronger with increasing continuum intensity. However, 
the anti-correlation between continuum intensity and tempera- 
ture in parts of the line forming layers leads to a smaller varia- 
tion of the line strength. 

The (3D) average and the LHD model indicated by sym- 
bols in Fig. |3] point towards a potential problem in our de- 



lable 3. Sulphur abundances in Frocyon. Col. (1) is the wavelength or the line; col. (2) is the bquivalent Width; col. (3) is the 
sulphur abundance, A(S), according to the CO^'BOLD 3D model; col. (4) is A(S) from the <3D) model; col. (5) is A(S) from the 
ATLAS ID model; col. (6) is A(S) from the LHD ID model; col. (7) and (8) are the A(S) from fitting using a CO^BOLD 3D 
and a ATLAS+SYNTHE grid, respectively; in all ID models a micro-turbulence of 2.1 kms ' was assumed; col. (9) is the 3D 
correction; col. (10) is the NLTE correction from Takeda et al. (2005a); col. (11) is the A(S) we adopted; col. (12) is the statistical 
error, q-mlt is of 1 .50 for the LHD model and 1 .25 for the ATLAS model. 
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Fig. 3. Procyon: joint probability density function (not nor- 
malised) of continuum intensity and (disk-centre) EW for the 
1045.9 nm line in the 3D model (grey contours). The labelled 
symbols mark the 3D and (3D) average, as well as the result 
for a LHD model with q'mlt=1-5. The micro-turbulent velocity 
for the ID models is 2.1 kms 



termination of 3D abundance corrections for Procyon. None 
of the about 6 x 10"* vertical profiles entering the calculation 
of the probability density corresponds to the (3D) model - at 
least as far as it concerns the line formation properties. The 
same holds for the LHD model. It appears plausible that the 
micro-turbulence prescribed in ID line formation calculations 
- motivated from observations - is too large in comparison to 
the eff'ective micro-turbulence intrinsic to the 3D model. We 
note he r e that our adopted micro-turbulence of 2. 1 km s"' from 
lSteffen(1985') is in agreement, within quo t ed errors, wit h those 
ofiuhrmann ('l99 8h;lGratton etal.l (ll996l);lTakeda et al.l(fl998 ) 
and .Allende Prieto et all (l2002h . while lTakeda et alj (Il996h pre- 
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Fig. 4. Sun: same as Fig. [3] for the solar case. The micro- 
turbulent velocity for the ID models is l.Okms 



fer a lower value of 1.4kms . Test calculations have shown 
that a reduction on ^niicro from 2.1 kms"' to 1.5 kms"' leads 
to a reduction of the abundance correction 3D-LHD by a fac- 
tor of two. In part, this may ex plain the diff'erent 3 D correc- 
tions we obtain in comparison to lNissen et al.l (120041) who find 
negligible corrections for an only slightly cooler atmosphere 
(reff=6191 K, log §=4.04, [M/H]=0.0). On the other hand, all 
sulphur lines which span a range of a factor five in equivalent 
width provide a consistent abundance in 3D when ID-NLTE 
effects are included. We cannot provide a resolution here, but 
the issue of the appropriate micro-turbulence in 3D- ID com- 
parisons clearly needs further investigation. 

In LTE we obtain a sulphur abundance of 7.400 + 0.131, 
while including NLTE corrections it becomes 7.278 + 0.029. 
Once NLTE corrections are applied the scatter becomes tiny 
and fully compatible with the noise in the data. We commented 



berore that the 869.3 nm hne is blended by molecular lines. Ir 
we discard this line A(S)=7.250 + 0.026. Thus, the difference 
is small. This is perhaps not surprising, since Procyon is hotter 
than the Sun, and molecular lines are less prominent in this 
spectral range. 

7.3. HD 33256 

The stellar parameters found in the literature for this star show 
little scatter We assume Teff = 6 454 K, which we derived by 
fitting the Ha wings with a grid of synthetic spectra computed 
with ATLAS +SYNTHE. We recall that the SYNTHE code 
computes the van der Wa als broadening of the Balmer lines 
according to the theory of Ali & GrierrJ ( 1965.) . and the Stark 
broadening according to lVidal et alJ (TqTsI). H ad we adopted 
the broadening theory of iBarklem et al. I (l2000l) we would have 
obtained a lower effec tive temperature, as pointed out e.g. by 
Bonifacio et al.l (l2007h . We assumed the stellar parameters log 
g = 4.00, VsinO> lOkms-i, ^^i^-o = 1.55kms-i. From 20 
Fe I lines we obtain [Fe/H] = 7.25 + 0.1 1; 18 lines of Fe ii give 
[Fe/H] = 7.29 ± 0.15. For the ranges at 600 nm and 800 nm 
we measure a S/N ratio of about 700, and S/N = 25 around 
1045 nm. 

The results of our analysis are listed in Table H] The errors 
of A(S), derived from these S/N ratios according to the formula 
(1), are given in the last column of the table. 

The LTE s ulphur abundance is 7.022 + 0.183; applying the 



For this star the contribution of the iron line to the 
1045.5 nm line of Mult. 3 is large. So considering only the 
869.4 nm and the 1045.6 nm line, the sulphur abundance is 
7.297 + 0.106, and 7.248 + 0.070 if we apply the NLTE 



iTakeda et alJ (12005a) NLTE corrections A(S)=6.900 + 0.091. 
After applying the NLTE corrections the line to line scatter is 
fully compatible with the expected measurement errors. The 
sulphur abundance with NLTE correction from the lines of 
Mult. 3 is 6.938 + 0.123, to be compared with A(S)=6.862 ± 
0.033 if we consider the 675.7 nm, 869.3 nm and 869.4 nm 
lines, or A(S)=6.861 ± 0.047 if we consider only the two lines 
of Mult. 6. As can be seen, the sulphur abundances obtained 
from the different lines are in agreement within one standard 
deviation. 



7.4. HD 25069 



We adopted the stellar parameter from lValenti & Fischeii (l2005h 
(T^eff = 4994K, logg = 3.53, and [Fe/H]=+0.10). 

The fit of the 675.7 nm triplet gives A(S)=7.249, but, due 
to the presence of a distortion in the line profile, the fit was 
performed in a range not including the blue wing of the line. 

S/N is 500 for the 600 nm and 800 nm ranges. In the range 
1045 nm we measure S/N=28. EWs and sulphur abundances 
are reported in Table |5] The errors according formula ( 1 ) are 
provided in the last column of Table |5] 

For this star the IR line at 1045.9 nm appears to have an 
unusual shape, the core appears broad and flat, unlike the other 
atomic lines in this spectrum, for this reason we rejected this 
line. The triplet at 675.7 nm should also be rejected because 
the line profile appears distorted, possibly by a cosmic ray hit. 



corrections of ITakeda et alJ (l2005al) . Considering all lines 
A(S)=7.254 + 0.102 in LTE, A(S) ^7.207 + 0.120 with NLTE 
corrections of' Takeda et al. (2005a). 

Since A(S)=7.199 from the 869.4 nm line, and 
A(S)=7.169 ± 0.182 from the the two lines considered in 
Mult. 3, we can conclude that the abundance determinations 
from the two lines are in good agreement. We attribute the 
abundance difference of 0.26 dex between the lines at 1045.6 
and 1045.9 nm to the low S/N of the spectrum of Mult. 3. 

7.5. € Eri (HD 22049) 

For this star the stellar parameters found in the litera- 
ture ar e again in good ag reement. We take the parame- 
ters of ISantos et"aD t004 (T^s = 5 073K, logg = 4.42, 
^micro=l-05kms~'and [ Fe/H] =-0.1 3), because these are also 
the parameters used by lEcuvillon et alJ (120041) for their sul- 
phur abundance determination. Th e projected rotational ve loc- 
ity y sin(/)=3.0km s ' stems from lNordstrom et alj ( 120041) . 

For this star we measure S/N=300 at 670 nm and a bit 
smaller than 300 at 870 nm, but S/N =50 at 1045 nm. The results 
are reported in Table |6] The 3D abundance correction given in 
column (5) is the one for a star with T^ff - 5000 K, logg=4.44 
and solar metallicity. 

Also for this star the contribution of the iron line to the 
1045.5 nm blend of Mult. 3 is high (12% of the whole line), 
therefore we exclude this line from the S abundance determi- 
nation. The error of A(S) according to formula (1) can be found 
in Table|6] last column. 

The LTE sulphur abundance is: A(S)=7.208 ± 0.048, with 
NLTE correction 7.185 + 0.067. A(S)=7. 148 ± 0.06 1 if we con- 
sider the two lines of Mult. 3, to be compared with 7.241 ±0.010 
from the 675.7 nm and 869. 4nm lines. For this star the agree- 
ment of the S abundance derived from the lines of Mult. 3 with 
that derived from other lines is worse than for other stars. We 
note also that the lines of Mult. 3 in the observed spectrum ap- 
pear somewhat distorted and a higher S/N spectrum would be 
desirable to verify this result. 

8. Discussion 



The main interest of this investigation is the level of concor- 
dance of the abundances derived from the lines of Mult. 3 with 
the abundances derived from lines of other multiplets. Our re- 
sults are summarised in Fig. |6] where the value of A(S) for 
each line is plotted for all stars. There is no evident trend of 
A(S) with respect to the line used to obtain the abundance, ex- 
cept for the Sun, for which we have noted the strong A(S)-EW 
correlation. At this stage it is not clear if accounting for NLTE 
effects in the 3D model may solve this problem. However, the 
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Table 4. Sulphur abundances in HD 33256. Col. (1) is the wavelength of the line; col. (2) is the Equivalent Width; col. (3) is A(S) 
from ATLAS+Linfor3D; col. (4) is A(S) from LHD+Lin for3D; col. (5) is A(S) from fitting using an ATLAS+SYNTHE grid; 
col. (6) is the NLTE correction from lTakeda et al. I (l2005ah : col. (7) is the adopted A(S); col. (8) is the statistical error of A(S). 
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good agreement obtained for the other stars is encouraging and 
suggests that S i Mult. 3 is indeed a valuable abundance indi- 
cator. The fact that these abundances are consistent for stars of 
different spectral types suggests that all the oscillator strength 
are on the same scale and there is no systematic difference be- 
tween the different multiplets. We repeat that the error of the 
log gf for all lines used here is of the order of 50 %, and better 



laboratory or theoretical oscillator strengths would be highly 
desirable. 

There is a tendency to underestimate the EWs when the ob- 
served spectrum is broadened by an instrumental profile. Tests 
on strong solar lines, for which we have a very high S/N spec- 
trum show a lower value for the measurements of the EW of 
more than 3 %, when broadened to a resolution of R~ 80 000. 



lable 5. Sulphur abundances in HU 25069. Ihe rormat is identical to lab.|4J 
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Table 6. Sulphur abundances in e Eri (HD 22049). Col. (1) is the wavelength of the line; col. (2) is the Equivalent Width; col. (3) 
is A(S) from ATLAS +Linfor3D; col. (4) is A(S) from LHD+Linfor3D; col. (5) is A(S) from fitting using an ATLAS+SYNTHE 
grid; col. (6) is the 3 correction as 3D-LH D for a model of Teff = 5000 K, logg = 4.44 and solar metallicity; col. (7) is the 
NLTE correction from lTakeda et al.l ( l2005al) ; col. (8) is the adopted A(S); col. (9) is the statistical error of A(S). 



Wave EW 
(nm) (pm) 



A(S) 



3D- ID 
FIT 



A(S) 



(1) 


(2) 


ATLAS 
(3) 


LHD 

(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


675.7 


0.62 


7.162 


7.234 


7.080 


-0.0360 




7.234 


0.025 


869.4 


1.02 


7.183 


7.252 


7.047 


-0.0304 


-0.004 


7.248 


0.023 


1045.5 


5.81 


7.092 


7.143 




0.0268 


-0.047 


7.096 


0.035 


1045.6 


2.52 


7.175 


7.240 




-0.0056 


-0.025 


7.215 


0.082 


1045.9 


4.36 


7.111 


7.169 


7.047 


0.0033 


-0.037 


7.132 


0.057 



When dealing with an observed spectrum of lower quality we 
expect this effect to be even larger. According to our simula- 
tions an IR low quality (S/N=20) spectrum can lead to measure 
an EW of a strong line which is 10% smaller than the true 
value, because the wings are lost in the noise. 

Our main conclusion is that Mult. 3 can be used success- 
fully to measure the sulphur abundance in the sample of stars 
we have considered in this work. We believe that this analysis 
can be extended to metal-poor stars, where the measurement 
of Mult. 3 lines with the IR spectrograph CRIRES seems 
particularly promising. 
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Appendix A: Remarks on individual stars 

1. Sun: 

The 675.7 nm sulphur triplet is not blended and it is well 
reproduced in comparison with synthetic spectra, when us- 
ing the atomic data reported in Table [1] We do not con- 
sider the sulphur triplet at 674.3 nm whose shape in the so- 
lar observed spectra is not well reproduced, perhaps due to 
blending by a CN line. The 674.8 nm sulphur triplet is also 
discarded because of a blend with vanadium, calcium, iron 
and titanium whose atomic data are not well known. 
The two sulphur lines of Mult. 6 in the range 870 nm are 
both well reproduced by synthetic spectra. The lines are 
situated very close to each other, so in the ID analysis we 
fitted both simultaneously. In the 3D analysis we concen- 
trated on the 869.4 nm line. The 869.3 nm line is blended 
with molecules (CN and C2), and it is weaker than the 
879.4 nm line. We computed the EW of the contribution 
of molecules in the range of this line to derive the sulphur 
abundance from this line. 

We do not consider in this work the Mult. 1 (920 nm) sul- 
phur lines, because they are affected by telluric lines. The 
922.8 nm line is the cleanest among the three lines, but it 
lies in the wing of the Paschen ^ H-line, so the abundance 
analysis is not free from possible systematic errors. 

2. Procyon: 



Allend e Prieto et al. ( I2002f) and later lAufdenberg et al 



( 120051) investigated Procyon's limb darkening with 3D 
models. Both groups found that 3D models predict 
a smaller degree of limb darkening than ID models. 
Aufdenberg et al. showed that an approximate overshooting 
introduced in ID Phoenix or ATLAS models can largely 
eliminate these differences. The case of Procyon is fortu- 
nate because it is a well studied visual binary system. The 
mass of the pri marv is well determined (see Gatewood & 
Han 2006 and Girard et al.ll2000l for the most recent astro- 
metric studies of this star). 

The angular diameter has been m easured directlv bv Brown 
et al. ( 1967i); .Hanburv Brown et al. ( 19 7Ji) and, more re - 
centlv. bv lShao et all (ll988h : lMozurkewich et all (Il99lh : 



Idi Be nedetto ( 1998); Nordgren et alj fl2UUll): Mozurkewich 
et al. C2003 ); Kervella et al. (20040. 

The distan ce from Hipparcos and orb ital motion studies are 
reported in lGatewood &Hanl(l2006l) . 
The surface gravity of this star is well known from the or- 
bital data and diameter measurement s: loR^= 4.0 with an 



No detailed abundance analysis has been done tor this star, 
but many measure s of its radial velocity exist up to 2005. 



error l ess than 0.1, in agreement with l Allende Prieto et al 
( l2002h . who obtained log g= 3.96 + 0.02. 



The effective temperature depends on the method used to 
derive it. Critical reviews of values obtained by various 
authors with different methods are summarised in .Stefferh 
( 1985 ) and in Kato et al. ( 1996h . The Teg values range from 
6 400 K (the lowest value derived from the continuum and 
IR fluxes) to more than 6 800 K (from the ionisation balance 
of Fe). By excluding the high values derived from line spec- 
trum analysis, the highest T^ff value depends on the value of 
the mixing-length adopted in the models; a lower mixing- 
length corresponds to a lower convective flux so to a higher 
temperature gradient at the bottom of the the atmosphere 
llCato & Sadakane 1982). A synthetic spectrum at an effec- 
tive temperature of 6500K with gMTT - 0.5 agrees with 
Balmer line profiles (IFuhrmann et al.ll 19941) . 
We give more weight to the value derived from the flux dis- 
tribution, adopt 6 500 K and ascribe the higher values re- 
quired to fit the ionisation balance to the inaccuracy of the 
structure of the atmosphere adopted in ID models. These 
models require a higher Tetf to describe the outer atmo- 
spheric layers where the lines are formed. 
A further proof of the inadequacy of the structure of the 
ID models is given by the extensive dis cussion on the 
derivation of the atmospheric parameters by iLuck & Heiter 
62005) based on the choice of the model that better fits the 
spectroscopic data; they obtain Teff - 6850K, logg = 4.55, 

fmicro — 2.4 km S 

A selection of stellar parameters from the hterature is given 
in Table lO 

3. HD 33256: 

HD 33256 is slightly cooler than Procyon (F5 IV-V) in spite 
of the earlier spectral type given in Simbad and the Bright 
Star Catalog (BSC), F2V. It is sUghtly metal deficient and 
so has weaker lines than a solar abundance star This is 
at the origin of its earlier than Procy qn's spectra l type. In 
fact the new accurate classification bv lGrav et al.l (l2003h is 
F5.5V (kF4, mF2). 

It is not far from the Galactic plane, in fact its Galactic coor- 
dinates are: 208.83, -24.83. It is a thin disk star according 
to Bensbv et al. ( 2003). Th e S ab undance has been mea- 
sured bv lTakada-Hidai et al. ( 20021). T hese authors adopted 
[Fe/H] from lEdvardsson et al. 1 11993'), while they derived 
temperature and gravity (see Table IA!2T i. The LTE derived 
S abundance, based on the 8 693 nm and 8 694 nm lines, is 
7.29. A selection of stellar parameters is listed in Table lA.2l 

4. HD 25069: 



Valenti & Fischer! (l2005h made an extensive study of the 



spectroscopic properties of cool stars. From one Keck spec- 
trum they derived (with Kurucz models): Tgif - 4994K, 
logg = 3.53, [M/H] = 0.10, ysin(0= 3.3 kms"', RV 
= 39.2 kms"'. Abundances of Na, Si, Ti, Fe, Ni are also 
given. The star seems to be cooler than spectral type G9, as 
foun d in Simbad data base. W e keep the steUar parameters 



from Valenti & Fisched (120051) . 

It is remarkable that, although for such a cool star the 
Ho- profile is very little sensitive to effective temperature, 
from the fitting of the wings of Ha we obtain: T^s = 
4 726 K, which is in very good agreement with the temper- 
ature derived from the B-V colour: B-V=1.00 implies Tgif 
= 4700K. 
e Eri: 

e Eri (HD 22049) is a much studied star (795 papers since 
1983). It is a variable star of BY Dra type (BY Dra stars 
are flare stars with cool spots, which cause photometric 
variations during the rota tion of the star). It has two sus- 
pecte d planets: e Er i b ("Campbell et al.l 1988 ). e Eri c 
(Guillen & Thorndike 2002). Ben edict et al.l (l2006l) deter- 
mined the mass of e Eri b from HST and ground-based 
astrometric and RV data, modelled its orbit, confirmed 
the existence of this companion an d discuss the existence 
of the pos sible tertiary invoked by lOuillen & Thorndike 
(120021) and lOzernov et al. I (l2000l ). The star has a high level 
of magnetic activity inferred from chromospheric activ- 
ity consistent with a relatively young age, less than 1 Gyr. 
Observational and theoretical searches for the signature 
of planetary/brown dwarf objects in the s tructure of the 



dust d isk around this star are underway by [Benedict et al 



(l2006h . A selection of stellar parameters is listed in Table 
A3] 
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iable A.l. Frocyon: stellar parameters with rererence. 



T ft 

eft 

K 






km s"' 


S micTo 

km s"' 


V sin(0 


R pfprPTipp 












2.8 


Grav (1981j 


6500 ± 


80 


4.0 ±0.1 




2.1 ±0.3 


< 4.5 


Stelfen(1985) 


6605 




4.13 


-0.06 


2.23 




Gratton et ak (1996) 


6500 




4.00 




1.4 


3.3 


Takeda et al. (1996) 


6470 




4.01 ±0.10 


-0.01 ±0.07 


1.91 ±0.20 


2.6 ± 1.0 


Fuhrmann (1998) 


6640 




4.13 




1.8 


6.6 


Takeda et al. (19^ 


6530 ± 


50 


3.96 ± 0.02 




2.2 


3.16 ± 0.5 


Allende Prieto et al. (2002) 


6543 ± 


84 


3.975 ±0.013 








Aufdenbere et al. (2005) 



Table A.l. SteUar parameters of HD 33256 





logg 


[Fe/H] 




V sin(0 


Ref 


K 






km s ' 


km s"' 




6550 


4.09 


-0.34 


1.4 




Nissen 1981 


6270 


4.0 








Cleee etal. 1981 


6270 


4.0 


-0.26 


1.0 




Tomkin et al. 1985 


6440 


4.05 


-0.30 






Lambert et al. 1991 


6300 


3.60 


-0.45 


1.40 




Zhao & Magain 1991 


6400 


3.95 


-0.336 







£riel & Boesgaard 1992 


6442 


4.05 


-0.30 






^dvardsson et al. 1993 


6386 


4.10 


-0.30 


2.10 




King & Boesgaard 1995 


6385 


4.10 


-0.30 






Chen et al. 2001 


6440 


3.99 


-0.30 


2.3 




Takada-Hidai et al. 2002 


6411 


3.87 


-0.30 


1.5 




Grav et al. 2003 


6427 


4.04 


-0.30 


1.90 




Bensbv et al. 2003 
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Table A.3. e Eri (HD 22049): stellar parameters with reference. 





log,? 


[Fe/H] 




Vsin(i) 


Reference 


K 






km s"' 


km s"' 




5180 


4.75 


-0.09 


1.25 




Drake & Smith (1993) 


5 076 


4.50 


0.05 






Alonso etal. (19^ 




4.84 






2.1 


Allende Prieto et al. (2000) 


5 135 


4.70 


-0.07 


1.14 




Santos et al. (2001) 


5 117 








3 


Nordstrom et ak (2004) 


5 073 


4.43 


-0.13 


1.05 




Santos et al. (2004) 


4992 










Ramirez & Melendez (2005) 


5 177 


4.72 


0.06 


0.62 




Takeda et al. (2005b) 


5 200 


4.50 


-0.04 


0.70 




Luck & Heiter (2005) 



